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I.  INTRODUCTION 


A  submarine  can  cause  a  local  distortion  of  the  earth’s  naturally  occurring  magnetic 
field  and  thereby  produce  a  magnetic  anomaly.  Submarine-generated  magnetic  effects  arise 
from  the  permanent  and  induced  magnetization  of  the  ferromagnetic  material  constituting  the 
structure  and  machinery  of  the  submarine  and  from  galvanic  corrosion  currents  flowing 
through  the  submarine  and  the  surrounding  water. 

Permanent  magnetization  results  from  the  long-term  alignment  of  the  submarine  with 
the  earth’s  magnetic  field  and  is  enhanced  if  the  hull  is  stressed  through  vibration,  mechanical 
shock  or  submergence  and  subsequent  relaxation.  Over  short  time  periods,  the  permanent 
magnetization  is  independent  of  the  submarine’s  heading. 

Even  if  the  submarine  is  not  permanently  magnetized,  it  can  still  produce  a  local 
distortion  of  the  fairly  uniform  magnetic  field  of  the  earth.  This  distortion  can  be  described  as 
if  it  originates  from  a  temporary,  orientation-dependent,  induced  magnetization  of  the 
submarine. 

Electric  currents  produce  magnetic  fields.  Galvanic  corrosion  currents  arise  from 
differences  in  electric  potential  of  interconnected  dissimilar  metals  immersed  in  an  electrolytic 
solution  such  as  a  bronze  propeller  in  electrical  contact  with  a  steel  hull  protected  against 
corrosion  by  attached  sacrificial  zinc  blocks,  all  immersed  in  sea  water.  The  potential 
difference  between  the  propeller  and  the  zincs  is  about  1  volt  and  the  resulting  current  may 
range  from  1  to  100  amperes.  Because  of  variations  in  resistance  (as  a  function  of  propeller 
shaft  angle)  between  the  rotating  propeller  shaft  and  its  bearings,  this  current  may  be 
modulated  at  the  shaft  rate.  Thus  a  submarine  may  be  surrounded  by  steady  and  alternating 
fields  resulting  from  corrosion  currents. 

Each  of  the  foregoing  contributes  to  the  submarine’s  magnetic  signature,  which 
represents  a  potentially  rich  source  of  ASW-relevant  information.  An  understanding  of  such 
signatures  should  contribute  to  the  optimum  processing  of  signals  which  may  yield  improved 
target  acquisition  range,  improved  search  strategy,  target  classification,  and  heading  and  depth 
information. 

This  present  report  is  a  continuation  of  work  described  in  reference  (a).  In  reference 
(a),  equations  were  developed  for  describing  the  magnetic  anomaly  of  a  submarine  as  it  would 
appear  to  an  airborne  scalar  total-field  magnetometer  passing  through  it  at  constant  velocity. 

In  this  earlier  study,  the  submarine  was  treated  as  both  a  point  ferromagnetic  dipole  and  as  an 
extended  ferromagnetic  dipole  possessing  both  induced  and  permanent  magnetic  moments.  In 
addition,  expressions  for  the  magnetic  anomaly  from  the  steady  galvanic  corrosion  current 
(but  not  from  the  alternating  component)  and  the  combination  of  all  of  the  above  were  derived 
as  observable  by  an  airborne  sensor.  In  the  earlier  study  it  became  apparent  that  a  great 
portion  of  the  variability  of  the  signatures  stemmed  from  not  the  target-centered  anomaly  itself 
but  rather  from  the  path  that  the  sensor  took  through  the  anomaly.  Therefore,  it  was  decided 
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in  the  present  study  to  calculate  the  anomaly  in  target-fixed  coordinates.  This  approach 
revealed  the  lobe  structure  of  the  anomaly  which,  in  turn,  enabled  a  clearer  understanding  of 
the  signatures  that  an  airborne  magnetometer  produces.  In  addition,  in  the  present  study,  the 
temporal  and  spatial  dependencies  of  the  magnetic  anomaly  from  modulated  (alternating) 
galvanic  currents  were  determined  and  are  presented  here  in  both  target-fixed  and  sensor-fixed 
coordinates. 
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II.  FERROMAGNETIC  ANOMALY  OF  A  SUBMARINE 


Ferromagnetic  Moment  of  a  Submarine 

The  magnetic  source  strength  of  a  submarine  is  described  in  terms  of  its  magnetic 
moment,  which  is  a  vector  quantity  whose  magnitude  and  direction  depend  principally  upon 
the  magnitude  and  direction  of  the  local  magnetic  field  of  the  earth,  the  mass  and  shape  of 
ferromagnetic  material  constituting  the  submarine,  and  the  orientation  of  the  submarine  in  the 
earth's  field. 

In  reference  (b),  Peizer  has  provided  simple  formulas  for  calculating  the  induced  and 
equilibrium  magnetic  dipole  moments  based  on  the  shape  and  volume  enclosed  within  a  steel 
prolate  spheroid  approximation  of  a  submarine.  These  formulas  were  generalized  somewhat 
and  recast  in  reference  (a)  as  follows: 


Mve  =  -V  Hv  T  (1) 

Ml  =  Mu  +  Mlp  =  L  Hl  T  +  Mlp  =  L  Hh  T  cos(8— ct)  +Mlp  (2) 

Ma  =Mai+Map  =  AHaT  +  Map  =  A  Hh  T  sin(P-a)  +  Map  (3) 

in  which 

Mve  =  equilibrium  vertical  ferromagnetic  moment  (y  •  ft3) 

Ml  =  longitudinal  component  of  ferromagnetic  moment  (y  •  ft3) 

Mu  =  induced  longitudinal  moment  (y  ■  ft3) 

Mlp  =  permanent  longitudinal  moment  (y  •  ft3) 

Ma  =  athwartship  component  of  ferromagnetic  moment  (y  •  ft3) 

Mai  =  induced  athwartship  moment  (y  •  ft3) 

Map  =  permanent  athwartship  moment  (y  •  ft3) 

V  =  coefficient  for  equilibrium  vertical  component  of  moment  (ft3/ton) 

L  =  coefficient  for  induced  longitudinal  component  of  moment  (ft3/ton) 

A  =  coefficient  for  induced  athwartship  component  of  moment  (ft3/ton) 

Hv  =  vertical  component  of  the  earth's  field  (y) 

Hl  =  longitudinal  component  of  the  earth's  field  (y) 

Ha  =  athwartship  component  of  the  earth's  field  (y) 

Hh  =  horizontal  component  of  the  earth's  field  (y) 

T  =  submerged  displacement  of  the  submarine  (tons) 
a  =  submarine  heading  (degrees  true) 

P  =  local  geomagnetic  declination  (degrees;  east,  positive). 

Figure  1  illustrates  the  magnetization  of  a  submarine  in  the  earth's  magnetic  field  and 
the  corresponding  vector  components  of  the  submarine's  magnetic  moment. 


3 


Figure  1.  Magnetization  of  a  submarine  in  the  earth's  magnetic  field. 


In  equations  (1),  (2),  and  (3)  the  positive  directions  of  Mve,  Ml,  and  MA  are  vertically 
upward,  stem-to-bow,  and  portside-to-starboard,  respectively.  The  positive  direction  of  Hv 
is,  consistent  with  convention,  vertically  downward.  The  magnitude  of  the  resultant  moment 
is 


M  =  (Mve2+  Ml2+  Ma2)172.  (4) 

Peizer  provides  values  for  the  coefficients  of  magnetic  moments  for  six  classes  of 
submarines.  He  gives  a  single  value  of  V  =  5.1  ft3/ton  for  all  six  classes  but  six  values  for  L 
ranging  from  10.6  to  16.6  ft3/ton  and  a  calculated  value  of  L  =  11.5  ft3/ton.  He  gives  three 
values  for  A  (2.9,  2.9,  and  3.6  ft3/ton)  and  a  suggested  value  of  A  =  3.0  ft3/ton. 


Magnetic  Induction  from  a  Point  Ferromagnetic  Dipole 

At  distances  that  are  large  in  comparison  with  the  length  of  a  submarine,  its  magnetic 
field  may  be  assumed  to  originate  from  a  point  magnetic  dipole  centered  at  the  submarine. 

The  vector  field  B  produced  by  a  point  magnetic  dipole  is  given  by 

B  =  -M/R3  +  (3  M-  R)  R/R5.  (5) 

In  equation  (5)  and  the  following  equations,  bold  type  indicates  vector  quantities.  R  is  the 
radius  vector  from  the  dipole  to  the  field  point  P  at  which  the  magnetic  field  is  to  be 
calculated,  R  is  the  magnitude  of  R,  B  is  the  magnetic  induction  vector  at  P,  and  M  is  the 
magnetic  moment  vector  of  the  dipole.  An  inspection  of  equation  (5)  shows  that,  at  a  given 
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distance  R,  the  magnitude  B  is  twice  as  great  along  the  dipole  axis  (for  which  MR  =  MR)  as 
along  the  perpendicular  bisector  of  the  dipole  (for  which  M-R  =  0). 

For  some  purposes,  it  is  useful  to  rewrite  equation  (5)  in  scalar  form  as 

Bd=  (3cos20+  1)wM/R3  (6) 

in  which  BD  is  the  magnitude  of  the  magnetic  induction  from  the  ferromagnetic  dipole  and  0  is 
the  angle  between  the  positive  dipole  axis  and  the  direction  to  P.  In  equation  (6)  it  is  more 
readily  apparent  that,  at  any  given  orientation,  the  magnetic  induction  is  directly  proportional 
to  the  magnetic  moment  and  it  varies  inversely  as  the  cube  of  the  distance  from  the  dipole.  It 
should  be  observed  that,  along  the  dipole  axis  (i.e.,  for  0  =  0°),  equation  (6)  reduces  to 
BD  =  2  M/R3  whereas,  perpendicular  to  the  dipole  axis  (i.e.,  for  0  =  90°),  BD  =  M/R3. 

If  0M  is  the  angle  of  the  dipole  axis  relative  to  a  coordinate  system  whose  origin  is  at 
the  dipole  and  0R  is  the  angle  of  a  line  from  the  origin  to  field  point  P,  then  0  =  0r-0m- 
If  equation  (6)  is  solved  for  R,  one  obtains  the  range  RB,  as  a  function  of  angle,  for  which  the 
magnetic  induction  from  the  ferromagnetic  dipole  has  the  magnitude  Bd,  i-e., 

Rb  =  {[3  cos2(0r-0m)  +  l]m  M/Bd}w  (7) 

In  figure  2,  equation  (7)  is  plotted  in  two  dimensions  for  M  =  1.55  x  109  y  •  ft3  and 
Bd  =  0-03  y.  This  figure  may  be  thought  of  as  a  vertical  slice  through  a  three-dimensional 
shape  that  would  result  from  rotating  the  curve  about  the  dipole  axis. 


Elevation  above  target  (ft) 


—  Magnetic  Induction  Earth's  mag  Held  -  Moment  vector  axis 


Figure  2.  Vertical  section  of  the  locus  of  points  for  which  magnetic  induction  from  a  point 
ferromagnetic  dipole  of  moment  1.55  x  109  y  •  ft3  has  a  constant  value  Bd  =  0.03  y. 
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A  more  complete  representation  of  the  magnetic  induction  from  the  dipole  must  also 
take  into  account  the  direction  0B  of  the  B  vector  at  each  point  P.  It  can  be  shown  that 

0B  =  0r  +  arctan[0.5  sin(0R-0M)  /  cos(0r-0m)].  (8) 

Figure  2  shows  36  vectors,  all  of  magnitude  BD  =  0-03  y,  at  orientation  angles  0B  given 
by  equation  (8).  The  array  of  parallel  dotted  lines  in  figure  2  represents  the  direction  of  the 
earth’s  magnetic  field. 


Magnetic  Anomaly  from  a  Point  Ferromagnetic  Dipole 

Scalar  magnetometers,  such  as  those  commonly  used  for  air  ASW,  measure  the 
magnitude  of  the  vector  sum  of  the  magnetic  induction  vectors  of  the  earth  and  of  the 
submarine.  The  magnetic  induction  of  the  earth  is  typically  thousands  of  times  greater  than  of 
the  submarine  at  ranges  of  interest;  thus  the  magnetic  signature  of  the  submarine  represents 
only  a  tiny  perturbation  in  the  observed  field.  The  large  background  field  of  the  earth  is 
usually  removed  by  scalar  subtraction  and/or  filtering;  the  residue  is  the  submarine’s  magnetic 
anomaly.  In  Appendix  A,  equation  (5)  is  expanded  and  equations  are  given  for  calculating  the 
magnetic  induction  components  and  the  magnetic  anomaly  in  the  earth’s  field  at  any  point  P. 

It  is  shown  in  Appendix  H  of  reference  (a)  that  the  magnetic  anomaly  of  a  submarine  is 
given,  to  an  excellent  approximation,  as  simply  the  magnitude  of  the  component  of  the 
submarine’s  magnetic  induction  vector  in  the  direction  of  the  earth’s  magnetic  induction 
vector.  If  the  vectors  drawn  in  figure  2  are  examined,  it  is  seen  that,  in  the  upper  right  and 
lower  left  portions  of  the  figure,  they  are  aligned  parallel  to  the  downward-directed  earth’s 
field  lines  and  therefore,  at  those  two  locations,  add  in  a  maximum  positive  sense  to  the 
earth’s  field.  On  the  other  hand,  at  two  locations  lying  along  the  perpendicular  to  that 
direction,  the  vectors  are  anti-parallel  to  the  direction  of  the  earth’s  field  and  therefore  add  to 
it  in  a  maximum  negative  sense.  Thus,  the  anomaly  exhibits  two  positive  and  two  negative 
maximum  absolute  values.  In  addition,  there  are  four  intermediate  angular  positions  at  which 
the  magnetic  induction  vectors  from  the  submarine  are  perpendicular  to  those  of  the  earth; 
nulls  (of  the  magnetic  anomaly)  occur  at  these  locations  because  the  components  of  these 
vectors  in  the  direction  of  the  earth’s  field  are  zero. 

Figure  3  shows,  in  submarine-fixed  coordinates,  a  vertical  slice  through  the  anomaly 
inscribed  within  the  magnetic  induction  curve.  One  could  obtain  the  anomaly  curve  of  figure 
3  by  calculating  and  plotting  the  projections  of  the  vectors  along  the  earth’s  field  lines  at  each 
point  of  the  curve  of  figure  2.  Figure  3  represents  the  locus  of  points  in  a  vertical  plane  for 
which  the  magnetic  anomaly  is  0.03  y.  The  two  larger  lobes  correspond  to  positive  portions 
of  the  anomaly  and  the  two  smaller  ones  to  negative  portions. 
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-  Magnetic  anomaly  Magnetic  induction  Earth’*  field  axis 
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Figure  3.  Vertical  section  of  magnetic  anomaly  from  a  point  ferromagnetic  dipole. 


The  need  for  vector  addition  can  be  avoided  by  use  of  the  following  scalar  approach. 
Let  0E  be  the  direction  of  the  earth’s  field  and,  as  before,  let  0M  be  the  direction  of  vector  M, 
and  0B  the  direction  of  vector  R.  The  magnetic  anomaly  T  as  a  function  of  angle  0R  can  be 
obtained  by  regarding  the  two  terms  on  the  right  side  of  equation  (5)  as  separate  vectors,  the 
first  in  the  direction  of  M  and  the  second  in  the  direction  of  R.  By  multiplying  the  cosines  of 
the  angles  between  these  vectors  and  the  earth’s  field  vectors  (i.e.  cos(0m-0e)  and  cos(0r-0e), 
respectively)  one  determines  their  components  in  the  direction  of  the  earth’s  magnetic  field. 
These  components  are  then  added  algebraically.  Thus 

T  =  -(M/R3)  cos(0m-0e)  +  3  M/R3  cos(0R-0M)  cos(0R-0E).  (9) 

It  follows  that  the  locus  of  ranges  Rr  as  a  function  of  0R  for  any  given  positive  or  negative 
value  of  T  is  given  by 


Rr3  =  [3  cos(0R-0M)  cos(0r-0e)  -  cos(0m-0e)]  M/T .  (10) 


The  factor  cos(0R-0M)  arises  from  the  scalar  product  MR  in  equation  (5). 
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Actually,  equation  (10)  was  used  to  plot  the  anomaly  curve  shown  in  figure  3.  [Note: 
It  was  necessary  to  insert  both  positive  and  negative  values  of  T  into  equation  (10)  to  obtain 
both  the  positive  and  negative  lobes  of  the  anomaly.  However,  only  the  non-negative  values 
of  R3  that  resulted  were  selected  and  their  cube  roots  taken  to  determine  R  as  a  function 
of  Or.] 


The  angular  positions  of  the  maximum  values  of  Rr  can  be  determined  from  equation 
(10)  by  taking  the  first  derivative  of  Rr3  with  respect  to  0R,  setting  it  equal  to  zero,  and 
solving  for  0R.  Thus 

cos(0R-0M)  [-sin(0R-0E)j  +  cos(0R-0E)  [-sin(0R-0M)j  =  0.  (11) 

Equation  (11)  can  be  rewritten  as  the  sine  of  the  sum  of  the  two  angles  (0r-0e)  and  (0r-0m)- 
Thus 


sin[(0R-0E)  +  (0r-0m)]  -  sin(2  0r  -  0E  -  0m)  =  0.  (12) 

Equation  (12)  is  satisfied  if  (2  0R  -  0E-  0m)  =  n  n  radians,  in  which  n  =  0,  1,  2,  3.  Therefore, 
maxima  occur  at  angles 


0R  =  (0E+  0m  +  n  ri)/2.  (13) 

For  n  =  0,  it  is  seen  that  0r  =  (0E  +  0m)/2.  That  is,  the  angle  at  which  the  first  maximum  of  the 
anomaly  occurs  is  simply  the  average  of  the  angles  of  the  earth’s  field  vector  and  the 
submarine’s  magnetic  moment  vector.  In  other  words,  the  axis  of  the  anomaly  bisects  the 
angle  between  these  two  vectors  as  shown  in  figure  3.  Other  maxima  occur  and  are  equally 
spaced  at  angular  separations  of  n/2  radians  or  90°. 

The  value  of  Rr3  at  the  first  maximum  can  be  obtained  by  inserting  0R  =  (0E  +  0M)/2 
into  equation  (10).  Thus 

Rr3a»xi  =  {3  cos2[(0e-0m)/2]  -  cos(0E-0M)}  M/T.  (14) 

Similarly,  the  value  of  R3  at  the  second  maximum  can  be  obtained  by  inserting 
0r  =  (0E  +  0m  +  ft)/2  into  equation  (13).  Thus 

Rr3m«2  =  {3  sin2[(0E-0M)/2]  +  cos(0E-0M)}  M/(-T).  (15) 

The  negative  sign  in  equation  (15)  implies  a  negative  value  must  be  inserted  for  T. 

By  symmetry,  the  third  and  fourth  maxima  are  equivalent  to  the  first  and  second 
maxima,  respectively. 
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If  the  angle  (0e-0m)  between  the  earth’s  magnetic  field  vector  and  the  submarine’s 
ferromagnetic  dipole  moment  vector  is  small,  equation  (14)  reduces  to  Rr3m«i  ~  2  M/T  and 
equation  (15)  reduces  to  Rr3nu*2  «  M/(-r).  Thus,  subject  to  this  assumption,  the  absolute 
values  of  the  anomaly  at  its  maxima  are  equal  to  the  values  of  Bd  at  the  same  angular 
positions. 

If  the  anomaly  curve  of  figure  3  were  to  be  rotated  about  the  indicated  axis  of  the 
anomaly,  a  three-dimensional  representation  of  the  ferromagnetic  anomaly  which  may 
surround  a  submarine  would  result  having  the  appearance  of  figure  4.1  It  consists  of  two 
parts:  a  doughnut-shaped  ring  or  toroid  and  what  appears  to  be  a  dumbbell  or  a  partially 
inflated  balloon  that  has  been  shoved  halfway  through  the  hole  of  the  toroid.  These  curved 
surfaces  are  the  locus  of  points  for  which  the  magnetic  induction  is  0.03  y  greater  than  or  less 
than  the  ambient  field  of  the  earth.  The  value  on  the  toroid  surface  is  negative  relative  to  the 
earth's  field  and  that  on  the  dumbbell  surface  is  positive. 


Elevation  above  target  (ft) 


4000  - 

3000 

2000  ! 

1000  - 

North 
o  - 

-  1000  ■- 

-2000  - 

-3000  !- 

-4000  - 

-6000  -4000 


-2000 


2000 


South, 


4000  6000 


Range  south  of  target  (ft) 


Figure  4.  Three-dimensional  plot  of  the  anomaly  from  a  point  ferromagnetic  dipole. 


Figure  5  illustrates  a  single  vertical  slice  through  the  anomaly  of  figure  4.  In  addition, 
it  shows  contours  for  +/-0.003  y,  +/-0.3  y,  and  +/-1.2  y. 


1  Actually,  the  data  used  in  preparing  figure  4  were  obtained  by  computing  the  magnetic  anomaly  as  it  would 
appear  to  a  scalar  magnetometer  transported  on  a  large  number  of  radial  approaches  toward  the  ferromagnetic 
dipole  and  determining  the  ranges  at  which  the  absolute  value  of  the  anomaly  equaled  0.03  y. 
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Elevation  above  target  (ft) 


Figure  5.  Vertical  section  of  a  ferromagnetic  anomaly  in  submarine-fixed  coordinates. 

Figures  6  and  7  illustrate  calculated  signatures  for  a  ferromagnetic  dipole  as  they  might 
appear  to  an  airborne  magnetometer  at  350-ft  altitude  passing  through  the  anomaly  shown  in 
figures  4  and  5.  Figure  6  shows  signatures  from  four  simulated  crossings  on  easterly  courses 
over  the  path  of  the  submarine  at  150-ft  depth  on  a  northerly  heading.  Curves  A  and  B 
correspond  to  passes  3000  ft  and  1500  ft  north  of  the  submarine  whereas  curves  C  and  D  are 
the  counterpart  signatures  south  of  the  submarine.  The  principal  feature  of  curve  A  is  a  weak 
positive  value  near  the  point  of  closest  approach  to  the  submarine  corresponding  to  passage 
near  the  upper  edge  of  the  forward,  downward-directed  lobe  of  the  anomaly.  The  path 
associated  with  curve  B  is  one  that  took  the  simulated  sensor  first  through  the  negative 
toroidal  portion  of  the  anomaly,  then  through  the  positive  dumbbell  portion,  and  then  out 
through  the  eastern  portion  of  the  negative  toroid.  Curves  C  and  D  exhibit  greater  amplitudes 
than  A  and  B  because  they  correspond  to  passages  near  the  center  of  the  upward-directed, 
positive  lobe  of  the  dumbbell.  Figure  7  shows  two  signatures;  the  larger  one  corresponds  to 
an  easterly  pass  directly  over  the  submarine  through  the  negative  toroidal  part  of  the  anomaly; 
the  other  represents  a  pass  on  a  northerly  course  1000  ft  to  the  west  of  the  submarine.  In  this 
latter  case,  the  sensor  passes  first  through  the  positive,  upward-directed  lobe  of  the  dumbbell 
and  then  through  the  negative  toroid. 
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Figure  6.  Calculated  MAD  signatures  for  sensor  aircraft  on  easterly  courses 
forward  and  aft  of  a  northbound  submarine. 


Magnetic  anomaly  (gamma) 
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Figure  7.  Calculated  MAD  signatures  for  sensor  aircraft  on  090°  and  000° 
passes  over/abeam  a  northbound  submarine. 


11 


Reference  (a)  presents  more  than  one  hundred  plots  of  calculated  submarine  anomalies 
as  they  would  appear  to  an  airborne  MAD  for  various  aircraft  and  submarine  headings  and 
encounter  geometries. 

Figure  8  is  a  computer-generated,  three-dimensional  anomaly  from  a  point 
ferromagnetic  anomaly  representing  the  locus  of  points  for  which  F  =  +/- 1  y  superimposed  to 
scale  on  a  drawing  of  a  submarine. 


Figure  8.  Three-dimensional  anomaly  surrounding  a  submarine. 

The  curved  surfaces  represent  the  locus  of  points  for  which 
T  =  +/-1  y. 


Figure  9  provides  three  views  of  a  computer-generated  plot  of  the  ferromagnetic 
anomaly  which  would  surround  a  submarine  on  a  north  magnetic  heading  in  the  Strait  of 
Hormuz  where  the  inclination  of  the  earth's  magnetic  field  is  40.3°  (down).  The  coordinate 
system  is  centered  at  a  point  ferromagnetic  dipole  that  is  assumed  to  be  the  magnetic 
representation  of  the  submarine.  The  submarine  was  assumed  to  possess  the  following 
components  of  magnetic  moment: 


Mve  =  -7.1  x  108  y  •  ft3 

(16) 

Ml  =  13.8  x  108  y  •  ft3 

(17) 

Ma  =0 

(18) 

M  =  15.5  x  108  y  •  ft3- 

(19) 
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The  depression  angle  of  the  resultant  moment  vector  can  be  obtained  from  equations 
(16)  and  (17)  as  arctan(MvE  /  ML)  =  212°.  Thus  the  submarine's  moment  vector  lies  at  an 
angle  that  is  between  the  submarine's  longitudinal  axis  (horizontal)  and  the  direction  of  the 
earth's  field  vector.  Figure  9(a)  provides  a  horizontal  view  looking  east;  figure  9(b)  is  a 
horizontal  view  viewing  north;  figure  9(c)  is  a  vertical  view  from  above. 

The  axis  of  highest  symmetry  of  the  anomaly  is  oriented  at  a  depression  angle  of 
33.75°,  which  is  mid-way  between  the  magnetic  moment  vector  of  the  submarine  (27.2°)  and 
the  earth's  field  vector  (40.3°).  A  magnetic  anomaly  detector  exhibiting  a  sensitivity  of  0.03  y 
would  provide  a  detectable  signal  anywhere  within  the  closed  volume  illustrated.  Note,  in 
figure  9(a),  that  a  pass  of  the  MAD  on  an  easterly  heading  directly  over  the  center  of  the 
submarine  would  fail  to  detect  it  if  the  vertical  separation  exceeds  about  2200  ft;  however, 
2000  ft  to  the  north  or  the  south,  detection  could  be  achieved  at  a  vertical  separation  of  3000 
ft.  On  the  other  hand,  it  would  be  possible  to  miss  a  detection  2500  ft  north  of  the  submarine 
at  a  vertical  separation  of  1000  ft  but  to  achieve  detection  by  flying  1000  ft  higher. 

Although  equation  (7)  describes  the  range  RB  at  which  the  magnetic  induction  from  a 
dipole  has  the  value  BD,  without  taking  into  account  the  earth's  field,  it  still  can  be  used  to 
estimate  the  radius  of  the  toroidal  portion  of  the  anomaly  by  setting  0r-0m  =  90°  and  to 
estimate  the  maximum  radial  extent  of  the  dumbbell  portion  by  setting  0R-0M  =  0°  and  solving 
for  RB.  Thus  a  good  approximation  of  the  radius  Rbt  of  the  toroid  is 

Rbt  =  -(M/Bd)173  =  -(15.5  x  108  y  ■  ft3  /  0.03  y)1/3  =  -3725  ft,  (20) 

and  a  good  approximation  of  the  maximum  radial  extent  RBd  of  the  dumbbell  is 

Rbd  =  (2  M/BD)1/3  =  (2  x  15.5  x  108  y  •  ft3  /  0.03  y),/3  =  4693  ft.  (21) 

The  negative  sign  in  equation  (20)  may  be  related  to  the  negative  sign  of  the  anomaly  within 
the  toroid. 

In  Appendix  B,  36  different  horizontal  views  of  the  anomaly  illustrated  in  figure  9  are 
shown  as  a  succession  of  stereo  pair  images  at  10-degree  increments  of  angle. 

In  reference  (a)  it  was  pointed  out  that  the  anomaly  from  a  ferromagnetic  dipole  as 
observed  with  a  scalar  magnetometer  reveals  little  information  about  the  heading  of  the  target 
submarine.  For  example,  it  was  shown  that  if  the  submarine  has  no  horizontal  permanent 
component  of  magnetic  moment,  the  maximum  horizontal  deviation  of  its  moment  from 
magnetic  north  would  be  about  36°.  Because  the  anomaly  axis  bisects  this  angle,  the 
maximum  deviation  of  the  anomaly  would  be  18°.  If  it  is  assumed  that  the  submarine  can 
possess  a  horizontal  component  of  permanent  magnetic  moment  that  does  not  exceed  about 
25%  of  the  maximum  induced  longitudinal  moment,  the  maximum  deviation  angle  as  a 
function  of  submarine  heading  will  not  exceed  50°  and  the  anomaly  will  not  deviate  by  more 
than  25°.  Thus,  it  can  be  concluded  that  the  anomaly  axis  of  a  submarine  with  no  active 
signature  modification  equipment  will  deviate  between  18°  and  25°  from  magnetic  north  as  the 
submarine  heading  changes  over  a  full  360°.  Thus,  the  ferromagnetic  anomaly  provides  little 
unambiguous  information  regarding  the  heading  of  the  submarine. 
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Figure  9.  Computer-generated  plots  of  a  ferromagnetic  anomaly  in  submarine-fixed 
coordinates. 

(a)  Horizontal  view  from  portside  viewing  east. 

(b)  Horizontal  view  from  astern  viewing  north. 

(c)  Vertical  view  from  above. 
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Magnetic  Induction  and  Anomaly  from  an  Extended  Ferromagnetic  Dipole 

For  observation  distances  that  are  not  large  in  comparison  with  the  length  of  a 
submarine,  the  point  dipole  assumption  breaks  down.  For  handling  such  cases,  it  was 
assumed  in  reference  (a)  that  a  better  representation  could  be  obtained  by  representing  the 
target  magnetically  as  comprising  two  magnetic  poles,  an  N-pole  and  an  S-pole,  separated  by 
a  distance  equal  to  about  90%  of  the  length  of  the  submarine.  In  Appendix  D  of  reference  (a), 
expressions  were  developed  for  calculating  the  magnetic  induction  from  an  extended 
ferromagnetic  dipole  of  arbitrary  length,  location,  and  orientation.  It  was  shown  that  if  the 
distance  from  the  extended  dipole  to  the  observation  point  exceeded  several  target  lengths,  the 
two  approaches  yielded  equivalent  results.  However,  it  was  also  shown  that,  at  short 
distances,  such  as  those  used  in  making  airborne  measurements  of  submarine  magnetic 
moments,  considerable  error  could  be  introduced  by  applying  the  point  dipole  approximation. 
For  example,  in  the  case  of  a  point  dipole,  the  separation  of  the  extrema  of  the  magnetic 
anomaly  signature  is  a  function  of  the  target-sensor  distance  at  the  point  of  closest  approach. 
(This  is  consistent  with  the  lobe  structure  of  the  anomaly  as  depicted  in  figures  3,  4,  and  5.) 

On  the  other  hand,  for  an  extended  dipole,  the  separation  of  extrema  is  a  function  of  not  only 
the  observation  distance,  but  also  the  length  of  the  extended  dipole.  As  a  consequence,  when 
signatures  are  analyzed  from  the  point  of  view  of  a  point  dipole,  the  inferred  range  to  the 
target  and  its  inferred  magnetic  moment  appear  to  be  greater  than  they  really  are. 

To  compare  the  anomalies  calculated  for  point  and  extended  ferromagetic  dipoles, 
figures  10  and  11  were  plotted  of  vertical  slices  of  the  magnetic  anomalies  in  target 
coordinates.  As  before,  three-dimensional  figures  could  be  generated  by  rotating  these  figures 
about  the  anomaly  axes.  Except  for  scale,  figure  10,  for  a  point  dipole,  is  similar  to  figure  5. 
Note  that  for  the  point  dipole,  the  transitions  between  the  positive  lobes  and  the  negative 
lobes  occur  along  axes  that  pass  through  at  the  dipole  itself  located  at  the  origin  and  that  the 
thickness  of  the  toroid  is  zero  at  the  origin.  On  the  other  hand,  figure  11  shows  the 
corresponding  anomaly  for  an  extended  dipole  of  215.55-ft  length.  In  this  case,  the  negative 
"toroid"  has  a  thickness  of  about  210  ft,  or  97%  of  the  dipole  length,  at  the  origin. 

Figures  12,  13,  and  14  show  the  same  data  as  figures  10  and  11,  but  plotted  in  a  way 
that  permits  easy  comparisons  between  the  point  and  extended  dipole  anomalies.  Note  that 
the  distance  scales  are  different  in  each  of  these  three  figures. 
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Figure  10.  Range  at  which  the  anomaly  from  a  point  ferromagnetic  dipole 
equals  1, 10,  and  100  gamma. 
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Figure  11.  Range  at  which  the  anomaly  from  a  215.55-ft  long  extended  ferromagnetic 
dipole  equals  1, 10,  and  100  gamma. 
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Figure  12.  Ranges  at  which  the  anomalies  from  point  and  extended 
ferromagnetic  dipoles  equal  1  gamma. 
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Figure  13.  Ranges  at  which  the  anomalies  from  point  and  extended 
ferromagnetic  dipoles  equal  10  gamma. 
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Figure  14.  Ranges  at  which  the  anomalies  from  point  and  extended 
ferromagnetic  dipoles  equal  100  gamma. 
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III.  MAGNETIC  ANOMALY  FROM  THE  STATIC  ELECTRIC 
CURRENT  MOMENT  OF  A  SUBMARINE 


As  indicated  in  Section  I,  a  submarine  immersed  in  seawater  behaves  like  a  short- 
circuited  galvanic  cell.  Corrosion  currents  arise  from  differences  in  electric  potentials  among 
the  interconnected  dissimilar  metals:  copper  (in  the  bronze  propeller),  iron  (in  the  steel  hull) 
and  zinc  (in  the  hull-protecting  sacrificial  electrodes).  The  standard  electrode  potentials  of 
copper,  iron  and  zinc  are,  respectively,  +0.34,  -0.44,  and  -0.76  volt.  Thus  potential 
differences  of  1.1  V  can  exist  between  copper  and  zinc,  0.78  V  between  copper  and  iron,  and 
0.32  V  between  iron  and  zinc.  Generally  the  largest  currents  flow  between  the  propeller  and 
the  nearest  zincs,  giving  rise  to  a  conventional  current  which  flows  from  the  propeller,  through 
the  drive  shaft,  bearings,  and  hull  and  then  through  the  water  back  to  the  propeller.  However, 
currents  may  also  flow  in  the  opposite  direction,  particularly  in  the  forward  part  of  the 
submarine,  between  the  iron  in  the  hull  and  the  zinc  blocks.  Thus  the  pattern  of  electric 
currents  may  be  very  complicated. 

In  the  present  analysis  it  is  assumed  that  the  net  conventional  current  in  the  submarine 
is  forward  from  the  stem  toward  the  bow  and  the  resulting  magnetic  effects  can  be  described 
as  if  they  originate  from  a  current  source  whose  dimensions  are  small  compared  to  the 
observation  distance. 

In  Appendix  C,  equations  from  reference  (c)  are  adapted  to  the  problem  of  calculating 
the  magnetic  anomaly  of  a  submarine  arising  from  galvanic  corrosion  currents.  In  reference 
(a),  equations  were  developed  and  exercised  for  calculating  the  anomaly  as  it  would  appear  to 
an  airborne  scalar  magnetometer  passing  in  the  vicinity  of  the  submarine.  In  the  present 
report,  the  anomaly  is  presented  in  submarine-fixed  coordinates  for  both  air  and  seawater 
surroundings,  analogous  to  the  three-dimensional  results  given  in  Section  II  for  a 
ferromagnetic  dipole.  The  anomaly  from  a  ferromagnetic  dipole  is  indifferent  to  whether  the 
observation  is  made  in  air  or  seawater.  This  is  not  the  case,  however,  for  the  anomaly  from  a 
horizontal  electric  current  dipole  (HED)  assumed  to  represent  the  submarine  in  this  context. 

Three  sets  of  equations  are  given  in  Appendix  C,  one  for  the  observation  point  P  in  air 
(i.e.,  altitude  z  >  0),  another  to  cover  the  subsurface  situation  (i.e.,  altitude  z  <  0),  and  a 
third  to  cover  the  special  simple  case  of  both  z  =  0  and  HED  depth  -h  =  0.  The  third  set  can 
be  obtained  from  either  the  first  or  second  set.  It  will  be  seen  that  the  HED  anomaly  is  much 
more  complex  than  that  of  the  ferromagnetic  dipole  in  that  there  is  no  axis  of  rotational 
symmetry  and  that  the  structure  is  a  strong  function  of  the  orientation  of  the  HED  in  the 
earth’s  magnetic  field. 

In  figures  15  through  19,  the  equations  of  Appendix  C  (Cases  I  and  II)  were  used  to 
calculate  the  locus  of  points  for  which  the  magnetic  anomaly  T  from  a  100-m  deep,  static. 
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1000-A  •  m  HED  equals  +/-0.01  y.  The  assumed  location  was  at  25°N,  58°E  in  the  Gulf  of 
Oman  where  the  components  of  the  earth’s  magnetic  field  are 

BEy  =  34271  y  (north  component) 

Be*  =  662  y  (east  component) 

Bez  =  -26582  y  (vertically  upward  component). 

Calculations  were  performed  for  two  HED  headings:  000°  magnetic  and  090°  magnetic.  For 
each  HED  orientation,  +/-0.01-y  contours  were  plotted:  first,  a  horizontal  slice  at  sea  level 
(z  =  0)  and  then  a  succession  of  nine  or  ten  east-west  vertical  slices  at  various  distances  y 
north  and  south  of  the  HED. 

Figure  15  shows  the  +/-0.01-y  contours  for  the  000°  magnetic  heading  of  the 
submarine.  Figure  15,  a  plan  view  at  the  sea  surface,  shows  four  lobes:  positive  ones  to  the 
southeast  and  northwest  and  negative  ones  to  the  southwest  and  northeast.  Figures  16(a) 
through  16(i)  show  the  vertical  slices  at  distances  2600,  2000, 1000,  and  0  m  south  of  the 
HED,  and  200,  400,  500, 1000,  and  1400  m  north  of  the  HED.  Roughly  equal  amounts  of  the 
anomaly  appear  above  and  below  the  air-water  interface.  It  is  interesting  to  note  that  an 
MAD  aircraft  flying  on  a  000°  magnetic  course  directly  over  the  submarine  on  this  000° 
magnetic  heading  would  fail  to  detect  its  HED  anomaly. 


Distance  north  of  HED  (m) 


Distance  east  of  HED  (m) 


Figure  15.  Sea-level  horizontal  section  of  magnetic  anomaly  from  a  static  HED  on  000°magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

HED:  depth,  100  m;  electric  current  moment,  1000  A- m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  16.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  000°  magnetic 
heading  in  submarine- fixed  coordinates  showing  +/-0.01-y  contours. 

(a)  2600  m  south  of  HED  (b)  2000  m  south  of  HED 

HED:  depth,  100  m;  electric  current  moment,  1000  A  •  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  16.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  000°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

(c)  1000  m  south  of  HED  (d)  0  m  south  of  HED 

HED:  depth,  100  m;  electric  ament  moment,  1000  A  •  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  16.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  000°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-Y  contours. 

(e)  200  m  north  of  HED  (f)  400  m  north  of  HED 

HED:  depth,  100  m;  electric  current  moment,  1000  A  ■  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  16.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  000°  magnetic 
heading  in  submarine- fixed  coordinates  showing  +/-0.01-y  contours. 

(g)  500  m  north  of  HED  (h)  1000  m  north  of  HED 

HED:  depth,  100  m;  electric  current  moment,  1000  A  •  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  16.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  000°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

(i)  1400  m  north  of  HED 

HED:  depth,  100  m;  orientation,  000°  magnetic;  electric  current  moment,  1000  A  ■  m. 
Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 


Figure  17  shows  the  calculated  anomaly  as  it  would  be  observed  with  an  airborne 
sensor  moving  from  west  to  east  at  an  altitude  of  100  m  at  a  distance  of  400  m  north  of  the 
HED;  it  should  be  compared  to  figure  16(f).  Figure  17  may  be  thought  of  as  one  of  many 
traces  used  to  determine  the  values  used  in  plotting  figure  16(f).  It  shows  scale  lines 
indicating  the  +/-0.01-y  anomaly  levels,  with  markers  spaced  at  distance  intervals  of  150  m. 
The  anomaly  curve  crosses  the  +0.01-y  level  at  distances  -444,  -42,  +794,  and  +1630  m  east 
of  the  HED;  similarly,  it  crosses  the  -0.01-y  level  at  distances  -1630,  -794,  +42,  and  +444  m 
east  of  the  HED.  These  same  values  could,  in  principle,  be  read  from  figure  16(f)  by  drawing 
a  horizontal  line  at  the  100-m  elevation  level  and  noting  its  eight  points  of  intersection  with 
the  curves. 
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- HED  anomaly  (gamma)  +  *0.01  gamma  •  -0.01  gamma 

Figure  17.  Magnetic  anomaly  along  west-to-east  flight  path  from  a  static  HED  in  seawater. 
Sensor:  altitude,  100  m;  course,  090°  magnetic  400  m  north  of  HED. 

HED:  depth,  100  m;  orientation,  000°  magnetic;  electric  current  moment,  1000  A -m. 
Location:  25°N,  58°E.  Earth's  magnetic  Geld  dip  angle:  37.8°. 


Figure  18  shows  the  +/-0.01-Y  contours  for  the  090°  magnetic  heading  of  the 
submarine.  Figure  18,  a  plan  view  at  the  sea  surface,  shows  three  lobes:  a  positive  one  to  the 
south  and  negative  ones  to  the  northwest  and  northeast.  Figures  19(a)  through  19(j)  show  the 
vertical  slices  at  distances  3500,  3000,  2000, 1000,  500,  250,  and  0  m  south  of  the  HED,  and 
1000,  1400,  and  1600  m  north  of  the  HED.  For  this  geometry,  the  positive  lobe  of  the 
anomaly  appears  mostly  below  sea  level  and  the  negative  portion  above  sea  level. 

The  parameter  values  chosen  for  exercising  this  model  (i.e.,  p  =  1000  A  •  m  and 
T  =  +/-0.01  y)  are  rather  optimistic  from  the  detection  point  of  view.  Because  the  static  HED 
anomaly  exhibits  a  generally  inverse  square  law  dependence  on  range,  reducing  the  ratio  p/T 
by  a  factor  of  100  (e.g.,  to  more  realistic  values  of  p  =  100  A  •  m  and  T  =  +/-0.1  y)  would 
reduce  detection  ranges  to  about  10%  of  those  implied  in  figures  15  and  18. 

The  static  HED  anomaly  offers  two  significant  advantages  provided  the  sensitivity  of 
airborne  magnetometers  can  be  improved  to  exploit  it:  (1)  Because  of  the  inverse  square  law 
fall-off  with  range  (as  opposed  to  the  inverse  cube  law  fall-offs  for  the  ferromagnetic  dipole 
anomaly  and  for  the  alternating  HED  anomaly),  greater  detection  ranges  are  potentially 
available.  (2)  Because  the  static  HED  anomaly  is  much  more  strongly  dependent  upon  the 
heading  of  the  submarine  than  is  the  ferromagnetic  dipole  anomaly,  the  former  provides  a 
much  greater  opportunity  for  obtaining  tactical  information  than  the  latter. 
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Distance  north  of  HED  (m) 


Figure  18.  Sea- level  horizontal  section  of  magnetic  anomaly  from  a  static  HED  on  090°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

HED:  depth,  100  m;  electric  current  moment,  1000  A  ■  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 


Distance  above  sea  level  (m) 


Distance  east  of  HED  (m) 


Figure  19.  East- west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  090°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

(a)  3500  m  south  of  HED 

HED:  depth,  100  m;  electric  ament  moment,  1000  A  •  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  19.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  090°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-Y  contours. 

(b)  3000  m  south  of  HED  (c)  2000  m  south  of  HED 

HED:  depth,  100  m;  electric  current  moment,  1000  A  -  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  19.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  090°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

(d)  1000  m  south  of  HED  (e)  500  m  south  of  HED 

HED:  depth,  100  m;  electric  current  moment,  1000  A  ■  m. 

Location:  25°N,  58°E.  Earth’s  magnetic  field  dip  angle:  37.8°. 
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Figure  19.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  090°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

(f)  250  m  south  of  HED  (g)  0  m  south  of  HED 

HED:  depth,  100  m;  electric  current  moment,  1000  A  ■  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 


30 


Distance  above  sea  level  (m) 


Distance  above  sea  level  (m) 


Distance  east  of  HED  (m) 


Figure  19.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  090°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

(h)  1000  m  north  of  HED  (i)  1400  m  north  of  HED 

HED:  depth,  100  m;  electric  ament  moment,  1000  A  •  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Distance  above  sea  level  (m) 


Figure  19.  East-west  vertical  section  of  magnetic  anomaly  from  a  static  HED  on  090°  magnetic 
heading  in  submarine-fixed  coordinates  showing  +/-0.01-y  contours. 

(j)  1600  m  north  of  HED 

HED:  depth,  100  m;  electric  current  moment,  1000  A  •  m. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 


Equations  (C-8)  through  (C-13)  can  be  combined  to  yield  an  equation  for  calculating 
the  magnetic  anomaly  T  in  a  horizontal,  sea-level  plane  from  a  static  HED  also  at  sea  level: 

T  =  [100  ps  /( p2  Be)]  [Bex  sin(4>  -  0)  -  BEy  cos(<j>  +  0)  +  Bez  sin  <j>].  (22) 

Note  that,  for  this  special  case,  at  any  given  angle,  T  varies  as  1/p  . 

Equation  (22)  can  be  solved  easily  for  p  and  plotted  as  a  function  of  compass  bearing 
to  field  point  P  for  any  given  value  of  I\  Figure  20  is  the  locus  of  points  for  which  the 
magnetic  anomaly  from  a  1000  A  ■  m-HED  equals  +/-0.01  gamma  for  a  variety  of  HED 
headings.  The  first  panel  in  figure  20A,  for  which  the  HED  heading  is  000°  magnetic,  is  the 
counterpart  of  figure  15.  Note  that  figure  15  consists  of  four  distinct  closed  curves  which  are 
not  conjoined  whereas  figure  20A  consists  of  similar  closed  curves  which  are  conjoined  at  the 
origin.  The  reason  for  this  difference  is  that  the  HED  was  at  a  depth  of  100  meters  for  the 
former  but  at  the  air-sea  interface  for  the  latter.  There  is  a  similar  relation  between  the  fourth 
panel  of  figure  20A  and  figure  18  for  a  heading  of  090°  magnetic.  Note  that  a  small  positive 
lobe  appears  at  the  surface  when  the  HED  depth  is  zero  but  not  when  it  is  100  m.  (This  lobe 
disappears  at  an  HED  depth  of  61  m.)  Note  also  in  figure  20  that  the  HED  anomaly  rotates 
180°  for  360°  of  change  in  submarine  heading  and  that  the  relative  sizes  of  the  lobes  vary  with 
rotation,  with  the  southernmost  lobe  being  the  largest. 
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Figure  20A.  Magnetic  anomaly  at  the  sea  surface,  represented  as  +/-0.01-y  contours  in 
submarine-fixed  coordinates,  from  a  1000-A-  m  static  HED  at  the  sea 
surface  for  six  headings  of  the  submarine. 

Submarine  headings:  000,  030,  060,  090, 120,  and  150  degrees  magnetic. 
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Figure  20B.  Magnetic  anomaly  at  the  sea  surface,  represented  as  +/-0.01-Y  contours  in 
submarine-fixed  coordinates,  from  a  1000-A-  m  static  HED  at  the  sea 
surface  for  six  headings  of  the  submarine. 

Submarine  headings:  180,  210, 240,  270, 300,  and  330  degrees  magnetic. 
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The  inverse  square  law  range  dependence  of  the  anomaly  from  a  sea-level,  static-HED 
anomaly  observed  at  sea  level  is  illustrated  in  figure  21.  For  this  set  of  curves,  three  ratios  of 
electric  current  moment  and  magnetometer  sensitivity  were  assumed:  105,  104,  and  103  A  •  m/y, 
corresponding  to  combinations  of  source-strength/sensitivity  ratios  as  great  as 
1000  A  ■  m  /  0.01  y  to  as  small  as  100  A  •  m  /  0.1  y.  The  linear  dimensions  of  successive  lobes 
are  larger  by  a  factor  of  101/2  (3.16).  In  a  realistic  situation,  the  smaller  lobes  would  probably 
be  swamped  by  the  ferromagnetic  anomaly,  which  follows  an  inverse  cube  law  dependence; 
however,  at  the  longer  ranges  consistent  with  larger  HED  moments  and  greater  equipment 
sensitivities,  the  HED  anomaly  would  dominate. 


Distance  north  of  HED  (m) 


. 10*5  A-m/gamma  - 10*4  A-m/gamma  ~ ~  10*3  A-m/gamma 

Figure  21.  Magnetic  anomaly  in  submarine-fixed  coordinates  for  three  ratios  of  static  HED 
electric  current  moment  and  magnetometer  sensitivity:  10s,  104,  and  103  A  •  m/y. 


The  dependence  of  the  appearance  of  the  static  HED  anomaly  upon  geographical 
location  is  illustrated  in  figures  22,  23,  and  24,  which  correspond,  respectively,  to  magnetic 
latitudes  of  000°  (magnetic  equator),  045°  (mid-latitude),  and  090°  (near  the  north  magnetic 
pole).  In  each  case,  the  heading  of  the  HED  is  000°  magnetic  and  it  is  positioned  at  the  sea 
surface;  the  anomaly  is  that  observed  at  sea  level.  The  anomaly  may  be  visualized  as  the 
surface  expression  of  a  tilting  three-dimensional  figure  whose  southern  portion  emerges  from 
the  sea  with  increasing  latitude  while  the  northern  portion  submerges.  Near  the  north 
magnetic  pole,  only  the  previously  southernmost  lobes  appear  at  the  surface.  Near  the 
magnetic  pole,  the  anomaly  rotates  about  its  vertical  axis  at  the  same  rate  as  the  HED  heading 
changes  whereas  near  the  magnetic  equator,  it  rotates  at  half  that  rate. 
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Figure  22.  Horizontal  sea-level  slice  of  magnetic  anomaly  in  submarine-fixed  coordinates  from 
a  sea-surface  static  HED  for  earth’s  magnetic  dip  angle  of  0°  (magnetic  equator). 
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Figure  23.  Horizontal  sea-level  slice  of  magnetic  anomaly  in  submarine-fixed  coordinates  from 
a  sea-surface  static  HED  for  earth’s  magnetic  dip  angle  of  45°  (mid-latitude). 
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Figure  24.  Horizontal  sea-level  slice  of  magnetic  anomaly  in  submarine-fixed  coordinates  from 
a  sea-surface  static  HED  for  earth’s  magnetic  dip  angle  of  90°  (north  magnetic  pole). 
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IV.  MAGNETIC  ANOMALY  FROM  THE  ALTERNATING  ELECTRIC 
CURRENT  MOMENT  OF  A  SUBMARINE 


General 

As  indicated  in  Sections  I  and  III,  a  submarine  immersed  in  seawater  behaves  like  a 
short-circuited  galvanic  cell  which,  because  of  variations  in  resistance  (as  a  function  of 
propeller  shaft  angle)  between  the  rotating  propeller  shaft  and  its  bearings,  may  produce  an 
electric  current  that  is  modulated  at  the  shaft  rate.  Equations  for  calculating  the  magnetic  field 
intensity  from  a  harmonic  electric  current  dipole  immersed  in  a  semi-infinite  conducting 
medium  were  taken  from  reference  (c)  and  applied  to  the  situation  at  hand  in  Appendix  D. 

If  equations  (D-10)  through  (D-15)  for  the  alternating  HED  are  compared  to  the 
corresponding  equations  (C-2)  through  (C-7)  for  the  static  HED,  certain  similarities  and 
differences  are  observed.  For  example,  an  exponential  attenuation  factor  appears  in  the 
equations  for  the  alternating  HED  but  not  for  the  static  HED.  Also,  the  equations  for  the 
static  HED  tend  to  exhibit  an  inverse  square  law  dependence  of  magnetic  induction  upon 
range  whereas,  for  the  alternating  HED,  the  two  horizontal  components  of  magnetic  induction 
tend  to  show  an  inverse  cube  law  dependence  upon  range  while  the  vertical  component 
exhibits  a  generally  inverse  fourth  power  dependence  upon  range.  If  the  equations  for  the 
special  case  of  both  the  field  point  P  and  the  HED  being  at  the  air-water  interface  are 
compared  (i.e.,  (D-16)  through  (D-18)  for  the  alternating  HED  compared  with  (C-8)  through 
(C-10)  for  the  static  HED),  the  differences  in  range  dependence  can  be  seen  more  readily. 

Equations  (D-16)  through  (D-18)  can  be  combined  to  yield  an  equation  for  calculating 
the  magnetic  anomaly  T  in  a  horizontal,  sea-level  plane  from  an  alternating  HED  also  at  sea 
level: 


T  =  [100  pa  /  (yi  p3  Be)]  {BEx  [sin((j>  +  0)  +  sin(<J>  -  0)] 

+  BEy  [cos(<j>  -  0)  -  3  cos(<J>  +  0)]  +  BEz  [6  sin  <J>  /  (yi  p)]}.  (23) 

Note  that,  for  this  special  case,  at  any  given  angle,  T  varies  as  1/p3  for  the  contributions  from 
the  horizontal  components  and  as  1/p4  for  the  contribution  from  the  vertical  component. 
Accordingly,  equation  (23)  doesn’t  lend  itself  to  an  easy  explicit  solution  for  p  as  could  be 
done  with  equation  (22)  for  the  case  of  the  static  HED.  However,  equation  (23)  permits  the 
interesting  observation  that  near  the  magnetic  equator,  where  Bez  =  0,  T  varies  as  1/p3  while 
near  the  magnetic  polar  regions,  BEx  and  BEy  equal  zero,  and  therefore  T  varies  as  1/p4. 

ELFE  Profiles  in  Submarine-Fixed  Coordinates 

For  a  comparison  of  the  magnetic  anomalies  from  static  and  alternating  HEDs, 
equations  (22)  and  (23)  are  plotted  in  figure  25  on  the  same  scale  in  polar  form  as  a  function 
of  compass  angle.  The  curves  show  the  magnetic  anomalies  at  sea  level  at  a  range  of  1000  m 
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for  static  and  alternating  current  moments  ps  =  pa  =  1000  A  •  m.  Other  assumed  values  are 
frequency  f  =  1  Hz,  seawater  conductivity  a  =  4  S/m,  and  seawater  dielectric  constant  k  =  78, 
yielding  a  seawater  propagation  constant  yi  =  0.005619  m'1.  The  HEDs  were  assumed  to  be 
on  a  000°  magnetic  heading  at  the  sea  surface  at  geographical  coordinates  25°N,  58°E  where 
the  earth’s  dip  angle  is  37.8°.  Note  that  at  this  range,  the  alternating  anomaly  has  a  structure 
similar  to  the  static  anomaly  but  which,  at  this  range,  is  weaker  despite  the  very  large  value 
assumed  for  pa. 


Anomaly  (gamma) 


- Alternating  HED  - Static  HED 

Figure  25.  Cross  plots  of  magnetic  anomalies  vs  compass  angle  at  1000-m  range  at  sea 
level  for  1000- A  •  m  static  and  alternating  HEDs  at  the  air-sea  interface. 


The  orientation  of  an  alternating  HED  anomaly  varies  as  a  function  of  HED  heading  in 
a  way  that  is  similar  to  that  from  a  static  HED.  Figure  26  shows  three  polar  plots  of  the 
magnetic  anomaly  at  sea  level  as  a  function  of  compass  angle  at  a  constant  range  of  1000  m 
from  a  sea-level,  100-A  •  m,  alternating  HED  for  HED  orientations  of  000°,  045°  and  090° 
true.  A  frequency  of  4  Hz  was  assumed.  The  geographical  coordinates  were  taken  to  be 
25°N,  58°E.  The  dotted  and  dashed  circles  correspond  to  scale  values  of  0.001  and  0.002 
gamma,  respectively.  The  polarities  shown  for  the  lobes  correspond  to  the  situation  at  some 
instant  in  time  at  which  the  alternating  component  of  the  galvanic  corrosion  current  in  the 
submarine  is  directed  from  the  stem  toward  the  bow.  During  the  next  half  cycle,  the  polarities 
would  be  reversed.  If  the  sensor  were  to  move  rapidly  from  within,  say,  a  positive  lobe  to  a 
negative  lobe,  a  sudden  phase  shift  would  be  observed. 
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Figure  26.  Polar  plot  of  magnetic  anomaly  at  a  distance  of  1000  m  from  an  alternating  100-A  •  m 
HED  vs  compass  angle,  for  three  HED  orientations:  (a)  000°,  (b)  045°,  (c)  090°. 


Figures  27  and  28  illustrate  the  structure  in  space  of  the  magnetic  anomaly  from  an 
alternating  HED  analogous  to  that  shown  in  figures  18  and  19  for  a  static  HED.  Equations 
(D-10)  through  (D-15)  and  equations  (C-ll),  (C-12),  and  (C-13)  were  used  to  produce  these 
curves,  which  show  the  envelope,  in  HED  coordinates,  for  which  the  anomaly  equals  +/- 
0.001  gamma.  The  submarine  was  assumed  to  be  at  a  depth  of  30  m  on  a  090°  magnetic 
heading  at  geographical  coordinates  25°N,  58°E,  where  the  earth’s  magnetic  field  dip  angle  is 
37.8°.  Its  alternating  electric  current  moment  was  assumed  to  be  100  A  •  m  and  to  have  a 
frequency  of  4  Hz.  Figure  27  shows  the  locus  of  points  in  a  horizontal  plane  at  sea  level  for 
which  the  anomaly  is  +0.001  y  for  the  lobes  lying  along  the  north-south  axis  and  -0.001  y  for 
the  lobes  along  the  east-west  axis.  These  polarities  correspond  to  the  instants  in  time  at  which 
the  current  within  the  submarine  in  the  stern-to-bow  direction  is  greater  than  the  average  dc 
value;  these  polarities  reverse  during  each  following  half  cycle.  In  figure  27,  conditions  1  and 
2  given  in  Appendix  D  for  “Field  point  in  air”  are  satisfied  for  ranges  R»89  m;  thus  values 
near  the  origin  should  be  disregarded.  Figure  28  shows  vertical  profiles,  in  three  east-west 
planes,  (a)  800  m  south  of  the  HED,  (b)  in  the  same  plane  as  the  HED,  and  (c)  800  m  north  of 
the  HED.  The  solid  portions  of  the  curves  in  figure  28  represent  values  that  satisfy  all  of  the 
conditions  given  in  Appendix  D  for  both  the  above-water  and  below-water  cases,  whereas  the 
dotted  portions  represent  terra  incognita  for  which  the  conditions  are  not  satisfied 
resoundingly.  It  is  seen  that  the  underwater  values  appear  valid  only  at  the  shallow  depths 
farthest  away  from  the  HED.  Because  the  equations  used  to  calculate  the  underwater  values 
assume  that  the  energy  propagates  from  the  HED  up  to  the  surface,  over  a  horizontal  path, 
and  down  to  the  field  point  P,  figure  28  underestimates  the  underwater  ranges  in  the  vicinity 
of  the  HED.  Because  the  primary  interest  in  this  study  is  in  investigating  phenomena  that  will 
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yield  detection  ranges  for  airborne  sensors  that  are  longer  than  those  currently  achievable  with 
conventional  MAD,  these  shortcomings  are  not  considered  serious. 


-3000  -2500  -2000  -1600  -1000  -600  0  600  1000  1600  2000  2600  3000 


Distance  east  of  HED  (m) 

Figure  27.  Sea-level  horizontal  section  of  magnetic  anomaly  from  an  alternating  HED  on  090° 
magnetic  heading  in  submarine-fixed  coordinates  showing  +/-0.001-Y  contours. 
HED:  depth,  30  m;  electric  current  moment,  100  A  •  m;  frequency,  4  Hz. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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Figure  28.  East-west  vertical  section  of  magnetic  anomaly  from  an  alternating  HED  on  090° 
magnetic  heading  in  submarine-fixed  coordinates  showing  +/-0.001-Y  contours, 
(a)  800  m  south  of  HED 

HED:  depth,  30  m;  electric  current  moment,  100  A- m;  frequency,  4  Hz. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 


(a) 
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Figure  28.  East-west  vertical  section  of  magnetic  anomaly  from  an  alternating  HED  on  090° 
magnetic  heading  in  submarine-fixed  coordinates  showing  +/-0.001-y  contours, 
(b)  0  m  north  of  HED  (c)  800  m  north  of  HED 
HED:  depth,  30  m;  electric  current  moment,  100  A  •  m;  frequency,  4  Hz. 

Location:  25°N,  58°E.  Earth's  magnetic  field  dip  angle:  37.8°. 
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If  a  scalar,  total-field  magnetometer  of  adequate  sensitivity  and  response 
characteristics  were  to  be  flown  on  selected  straight  line  paths  through  the  alternating  HED 
anomaly  depicted  in  figures  27  and  28,  the  ELFE  signals  would  exhibit  envelopes  such  as 
shown  in  figures  29  and  30.  Figure  29  is  the  envelope  (or  the  motion-produced  modulation) 
of  the  ELFE  signal  resulting  from  the  movement  of  the  aircraft  past  the  HED  on  a  north 
magnetic  heading,  500  m  east  of  the  HED.  The  simulated  sensor  altitude  is  50  m.  Figure  29 
may  be  considered  as  the  amplitude  profile  of  the  anomaly  along  the  “500-m  east”  grid  line  of 
figure  27.  For  ease  in  making  the  comparison,  both  figures  have  the  same  distance  scales.  In 
figure  29,  two  horizontal  scale  lines,  calibrated  at  50-m  intervals,  indicate  where  the  curve 
crosses  the  +0.001-  and  -0.001-gamma  levels,  corresponding  to  the  six  intersections  of  the 
lobes  of  figure  27  with  the  “500-m  east”  grid  line.  (The  agreement  is  not  exact  because 
figure  27  is  a  sea-level  plot  whereas  the  altitude  for  figure  29  is  50  m.) 


Anomaly  (gamma) 


HED:  100  A-m,  30  m  deep,  090  deg,  4  Hz. 

Sensor:  50  m  alt,  000  deg  mag,  100  m/s, 

500  m  east  of  HED.  25N,  58E. 

Figure  29.  Envelope  of  magnetic  anomaly  in  air  from  an  alternating  HED  in  seawater 

corresponding  to  a  northerly  sensor  aircraft  pass  500  m  to  the  east  of  the  HED. 
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Figure  30  is  similar  to  figure  29  except  the  simulated  sensor  moves  on  an  easterly 
heading,  800  m  south  of  the  HED  and  matches  up  with  a  cut  through  figure  28(a). 


Anomaly  (gamma) 


HED:  100  A-m,  30  m  deep,  090  deg,  4  Hz. 
Sensor:  50  m  alt,  090  deg  mag,  100  m/s, 
800  m  south  of  HED.  25N,  58E. 


Figure  30.  Envelope  of  magnetic  anomaly  in  air  from  an  alternating  HED  in  seawater 

corresponding  to  an  easterly  sensor  aircraft  pass  800  m  to  the  south  of  the  HED. 
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Temporal  Characteristics  of  ELFE  Signatures 

The  ELFE  signature  appears  to  arise  from  abrupt,  periodic,  short-duration  increases  in 
electrical  resistance  between  the  propeller  drive  shaft  and  its  hull-mounted  bearings.  Even  if 
these  changes  in  resistance  are  assumed  to  be  instantaneous,  the  resulting  changes  in  electric 
current  will  be  smeared  out  in  time  because  of  the  inductive  and  resistive  properties  of  the 
system.  Because  the  ELFE  signals  are  “spikey,”  they  will  contain  many  frequencies  that  are 
harmonics  of  the  shaft-rate  fundamental.  As  the  ELFE  signal  propagates  through  sea  water, 
energy  at  the  various  frequencies  propagates  at  different  speeds,  contributing  to  a  further 
time-broadening  of  the  pulses.  In  addition,  attenuation  in  seawater  is  a  function  of  frequency, 
causing  further  changes  in  the  shape  of  received  pulses.  In  the  following  paragraphs,  an 
attempt  is  made  to  obtain  a  crude  estimate  of  the  shape  of  the  received  pulses. 

The  change  Ai  in  electric  current  resulting  from  a  small  change  AR  in  resistance  in  an 
electrical  circuit  of  emf  g,  resistance  R,  inductance  L  and  inductive  time  constant  Tl  =  L/R  as  a 
function  of  time  t  is 


Ai  =  -  g  (t  /  RL  +  1/R2)  exp(-R  t  /  L)  AR 


(24) 


or  approximately  as 


Ai  oc  -  g  exp(-t  /t  l)  AR.  (25) 

If  it  is  assumed  that  the  resistance  R  of  the  submarine-seawater  circuit  is  about  0.03  Q 
(consistent  with  an  emf  g  =  1  V  and  a  current  i  =  33  A)  and  the  inductance  L  =  1  microhenry 
(consistent  with  a  magnetic  induction  of  0.03  y  at  a  range  of  250  m  from  a  single-turn  coil 
carrying  a  steady  current  of  33  A)  the  inductive  time  constant  tl  =  L/R  =  10'6/0.03  =  33  ps, 
which  may  be  considered  the  intrinsic  rise  time  of  the  ELFE  pulse. 

The  speed  of  propagation  of  a  plane  electromagnetic  wave  of  frequency  f  in  seawater 
of  conductivity  4  S/m  is 


v  =  1600  f1/2  (m/s).  (26) 

If  the  depth  (-h)  of  the  alternating  HED  is  100  m  and  the  ELFE  fundamental  frequency  is  1 
Hz,  the  transit  time  to  the  surface  for  the  first  ten  harmonics  will  range  from  19.8  ms  to  62.5 
ms,  a  spread  of  42.7  ms. 

The  attenuation  coefficient  k  for  a  plane  electromagnetic  wave  in  seawater  of 
conductivity  4  S/m  is 

k  =  0.004  f1/2  (nT1)  (27) 
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and  the  transmittance 


t  =  exp(k  h)  =  exp(0.004  f1/2  h). 


(28) 


Thus,  if  h  =  -100  m,  transmittances  for  the  first  ten  harmonics  (i.e.,  for  1  Hz  <  f  ^  10  Hz,) 
will  range  from  0.67  to  0.28.  In  accordance  with  equations  (27)  and  (28),  the  lowest 
frequencies  will  suffer  the  least  attenuation  but,  in  accordance  with  equation  (26),  these 
lowest  frequencies  will  incur  the  longest  delays  in  propagating  through  the  seawater.  These 
lowest  frequencies  would  also  be  expected  to  have  the  highest  amplitudes. 

The  foregoing  considerations  were  applied  to  calculate  waveforms  for  ELFE  pulses 
after  they  emerge  from  the  sea  for  a  fundamental  frequency  of  1  Hz  and  HED  depths  of  1, 10, 
50,  and  100  m.  The  initial  waveform  was  assumed  to  have  a  very  short  duration  and  to  have  a 
sawtooth  shape  which  could  be  represented  by  a  Fourier  series  for  which  the  initial  relative 
amplitude  of  the  n'h  component  was  1/n.  Calculations  were  performed  out  to  the  thirtieth 
harmonic  and  the  results  are  shown  in  figure  31.  Widths  of  the  four  waveforms  at  their 
respective  half-amplitude  levels  are  0.18  ms,  1.83  ms,  8.6  ms,  and  16  ms  for  HED  depths  of 
1  m,  10  m,  50  m,  and  100  m. 


0.1  1  io  100 

Elapsed  time  (ms) 

— -  0  •  1  m  -4-  D  •  10  m  -%r  0  •  60  m  -B-  D  •  100  m 

Figure  31.  Calculated  waveforms  for  ELFE  pulses  for  alternating  HED 
depths  of  1  m,  10  m,  50  m,  and  100  m. 

If  a  sampling,  scalar,  total-field  magnetometer  is  used  to  observe  ELFE  signatures  it  is 
important  that  (1)  the  duration  of  each  sample  be  small  compared  to  the  duration  of  the  ELFE 
pulse,  and  (2)  an  adequate  number  of  samples  be  taken  during  each  ELFE  pulse.  If  the  time 
per  sample  is  too  long,  the  signal  will  be  diluted  out  over  time.  If  the  observation  time  per 
sample  is  appropriately  short  (i.e.,  much  shorter  than  the  ELFE  pulse  duration)  but  the 
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interval  between  samples  is  not  small  compared  to  the  ELFE  pulse  duration,  some  pulses  will 
be  missed  completely  or  only  a  small  fraction  of  the  individual  pulses  will  be  detected. 

The  AN/ASQ-208  takes  432  samples  per  second;  the  time  interval  between  successive 
samples  is,  therefore,  2.31  ms.  For  an  HED  depth  of  100  m,  the  received  ELFE  half¬ 
amplitude  pulse  duration  is  shown  in  figure  31  to  be  16  ms;  therefore,  on  average,  a  probably 
adequate  6.9  samples  would  be  taken  per  pulse.  For  the  shallower  depths,  on  the  other  hand, 
the  duration  of  the  pulses  may  be  shorter  than  the  time  interval  between  samples. 
Consequently,  many  of  these  short  duration  pulses  could  be  sampled  inadequately  or  even 
missed  completely.  It  appears  that  if  ELFE  pulses  from  alternating  HEDs  as  shallow  as  10  m 
are  to  be  sampled  at  a  rate  of  at  least  two  per  pulse,  the  interval  between  samples  should  be 
reduced  to  1.83/2  ms,  corresponding  to  1093  samples  per  second. 

Figure  32  illustrates  the  signature  that  an  ideal,  scalar,  total-field  magnetometer  would 
produce  upon  passing  on  an  easterly  heading  800  m  to  the  south  of  a  30-m  deep,  4-Hz, 
alternating  HED  on  a  090°  magnetic  heading.  (The  envelope  of  this  signature  appears  in 
figure  30  on  a  scale  of  distance  rather  than  time.) 


Anomaly  (gamma) 


HED:  100  A-m,  30  m  deep.  090  deg.  4  Hz. 

Sensor:  50  m  alt,  090  deg  mag,  100  m/s, 

800  m  south  of  HED.  25N,  58E 

Figure  32.  Magnetic  anomaly  in  air  from  an  alternating  HED  in  seawater  corresponding 
to  an  easterly  sensor  aircraft  pass  800  m  to  the  south  of  the  HED. 
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V.  COMBINATION  OF  FERROMAGNETIC,  STATIC  HED,  AND 
ALTERNATING  HED  ANOMALIES 


In  the  foregoing  sections,  each  of  the  types  of  anomalies  was  treated  separately. 
Because  each  anomaly  is  simply  the  projection  of  magnetic  induction  vectors  onto  the  earth’s 
magnetic  induction  vectors,  the  anomalies  at  any  point  can  be  combined  by  algebraic  addition. 
Figure  33  shows  anomaly  signatures  for  a  ferromagnetic  dipole,  a  static  HED,  an  alternating 
HED,  and  the  combination  of  the  three.  The  target  was  assumed  to  be  a  7600-ton  submarine 
operating  at  a  depth  of  30  m  on  a  270°  magnetic  heading  at  coordinates  25°N,  58°E.  Its  static 
electric  current  moment  was  assumed  to  be  1000  A  •  m,  directed  from  the  stem  toward  the 
bow.  The  alternating  electric  current  moment  was  taken  as  100  A  ■  m  from  stern  to  bow  and 
to  have  a  frequency  of  4  Hz.  The  sensor  aircraft  was  assumed  to  be  flying  at  a  speed  of 
100  m/s,  on  a  course  of  090°  magnetic,  at  an  altitude  of  50  m,  800  m  to  the  south  of  the 
target. 


Anomaly  (gamma) 


-  Static  HED  anomaly  -  Farromagnet  anomaly 

-  Combined  anomaly  Alternating  HED 

HED:  1000/100  A-m,  30  m  deep.  270  deg,  4 
Hz;  7600  tons;  Sensor-.  50  n  alt,  100  m/s 
090  deg,  800  m  south  of  HED,  25N,  58E 

Figure  33.  Anomaly  signatures  for  a  ferromagnetic  dipole,  a  static  HED, 
an  alternating  HED,  and  the  combination  of  the  three. 

See  text  for  assumptions. 


This  particular  combination  of  parameters  was  chosen  to  enable  each  of  the  effects  to 
be  seen  on  the  same  chart.  For  example,  the  point  of  closest  approach  to  the  target  was 
800  m,  which  is  beyond  typical  current  MAD  range.  At  significantly  shorter  ranges,  the 
ferromagnetic  dipole  anomaly  would  have  swamped  the  HED  anomalies.  The  westerly 
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heading  was  chosen  for  the  submarine  so  that  the  static  HED  anomaly  would  partially  cancel 
the  ferromagnetic  anomaly,  thus  increasing  the  relative  magnitude  of  the  alternating  HED 
anomaly.  For  the  assumed  conditions,  the  magnitude  of  the  static  HED  anomaly  exceeds  that 
of  the  ferromagnetic  anomaly.  This  occurs,  in  part,  because  the  ferromagnetic  anomaly 
decreases  as  the  cube  of  distance  whereas  the  static  HED  anomaly  decreases  as  the  square  of 
distance.  Obviously,  many  other  combinations  are  possible. 

The  particular  aircraft  path  assumed  for  figure  33  was  chosen  to  keep  the  frequency 
components  of  the  ferromagnetic  signal  low  in  comparison  to  the  alternating  HED  frequency 
to  enable  easy  separation  of  the  latter  from  the  former  by  simple  high-pass  filtering.  If,  for 
example,  the  aircraft  had  flown  at  low  altitude  directly  over  the  submarine,  with  both  on  000° 
headings,  the  higher  harmonics  of  the  ferromagnetic  signal  would  seriously  compete  with  the 
4-Hz  alternating  HED  signal  and  probably  render  it  undetectable. 

A  convenient  rule  of  thumb  for  estimating  the  fundamental  frequency  fo  of  the 
ferromagnetic  anomaly  can  be  inferred  by  referring  to  figure  5  and  noting  that  if  the  sensor 
were  to  move  in  a  circular  path  in  the  plane  of  the  figure  at  speed  v  at  a  distance  R  from  the 
origin,  it  would  sense  two  cycles  of  change  in  the  magnetic  field  per  360°  of  rotation.  Thus 

6>-v/(jiR)  (Hz).  (29) 

To  illustrate  equation  (29),  consider  a  case  in  which  the  sensor  aircraft  were  to  pass  on  a 
northerly  course  at  a  speed  of  100  m/s  and  an  altitude  of  50  m  over  a  30-m  deep  submarine. 
The  fundamental  frequency  t0  of  the  ferromagnetic  anomaly  directly  over  the  submarine  would 
be  0.4  Hz,  which,  with  its  harmonics,  approximates  the  ELFE  frequencies  associated  with  the 
propeller  rotation  rate.  The  magnitude  of  the  ferromagnetic  signal  would  be  many  tens  of 
gamma  and  thus  would  be  extremely  difficult,  if  not  impossible,  to  separate  from  the  ELFE 
signal. 
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VI.  DISCUSSION 


For  current  generation  submarines  and  current  generation  MAD  sensors  in  fixed-wing 
aircraft,  the  dominant  signature  arises  from  the  ferromagnetic  moment  of  the  submarine.  The 
anomaly  from  the  static  electric  current  moment  serves  to  introduce  small  distortions  of  the 
ferromagnetic  signature;  the  anomaly  from  the  alternating  component  of  the  electric  current 
moment  is  essentially  undetectable,  largely  because  of  competition  from  the  ferromagnetic 
anomaly.  Both  the  ferromagnetic  and  alternating  HED  anomalies  drop  off  largely  with  the 
cube  of  the  range,  whereas  the  anomaly  from  the  static  HED  drops  off  as  the  square  of  the 
range. 


If  MAD  sensors  of  significantly  improved  usable  sensitivity  can  be  developed  and 
competing  noise  can  be  driven  down,  detection  ranges  may  improve  to  the  point  where  the 
static  HED  anomaly  dominates.  At  these  longer  ranges,  the  ferromagnetic  anomaly  and  the 
alternating  HED  anomaly  will  have  dropped  off  in  the  same  proportion;  however,  because  the 
HED  frequency  is  independent  of  aircraft  location  but  the  frequency  of  the  ferromagnetic 
anomaly  varies  inversely  as  range,  it  may  be  posible  to  separate  the  two  signatures  at  the 
longer  ranges.  Detection  of  the  alternating  HED  signature  would  provide  information  on  the 
submarine’s  state  of  motion  that  the  other  signatures  do  not  provide. 

Use  of  the  static  HED  signature  (at  the  longer  ranges  achievable  at  higher  sensitivities) 
could  provide  information  on  the  heading  of  the  target  submarine  that  is  not  generally  (and 
unambiguously)  available  from  the  ferromagnetic  signature. 

Useful  detection  of  ELFE  signals  by  a  scalar,  total-field  magnetometer  in  a  fixed-wing 
aircraft  appears  improbable  as  compared  to  detection  of  a  passing  submarine  by  a  fixed-site 
sensor.  Fixed-site  sensors  offer  the  relative  advantages  of  freedom  from  platform  and 
geological  noise,  much  longer  integration  times,  the  ability  to  exploit  the  vector  nature  of  the 
signature,  and  a  much  larger  separation  in  frequency  between  the  ELFE  signal  and  the 
ferromagnetic  signal. 
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APPENDIX  A 


MAGNETIC  ANOMALY  FROM  A  POINT  FERROMAGNETIC  DIPOLE 

The  magnetic  induction  B  produced  at  a  field  point  P  by  a  point  ferromagnetic  dipole 
is  given  in  vector  notation  by  the  following  expression 

B  =  -  M/R3  +  3(M  •  R)  R/R5  (A-l) 

in  which  M  is  the  magnetic  moment  vector  of  the  dipole, 

R  is  a  range  vector  extending  from  the  dipole  to  point  P,  and 
R  is  the  magnitude  of  R. 

Let  P(x,y,z)  =  field  point  at  which  anomaly  T  is  to  be  calculated 
h  =  elevation  of  dipole  relative  to  sea  level  (-depth)  (ft) 
x  =  Rx  =  distance  of  P  east  of  dipole  (ft) 
y  =  Ry  =  distance  of  P  north  of  dipole  (ft) 
z  =  Rz  +  h  =  distance  of  P  above  sea  level  (ft) 

B*  =  easterly  component  of  magnetic  induction  from  dipole  (y) 

By  =  northerly  component  of  magnetic  induction  from  dipole  (y) 

Bz  =  vertically  upward  component  of  magnetic  induction  from  dipole  (y) 

Be*  =  easterly  component  of  earth’s  magnetic  induction  (nT  or  y) 

Bny  =  northerly  component  of  earth’s  magnetic  induction  (nT  or  y) 

Bez  =  vertically  upward  component  of  earth’s  magnetic  induction  (nT  or  y) 

Be  =  magnitude  of  the  earth’s  magnetic  induction  vector  at  P  (nT  or  y) 

R  =  [x2  +  y2  +  (z  -  h f]m 

Be  =  (Bex2  +  BEy2  +  Be*2)  1/2  (nT  or  y) 
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Each  of  the  vectors  in  equation  (A-l)  can  be  resolved  into  easterly  (x),  northerly  (y), 
and  vertically  upward  (z)  components  as  follows: 

B  =  i  Bx  +  j  By  +  k  Bz  (A-2) 

M  =  i  Mx  +  j  My  +  k  Mz  (A-3) 

R  =  i  Rx  +  j  Ry  +  kRz  =  ix+jy  +  k(z-h)  (A-4) 

Substituting  equations  (A-2),  (A-3),  and  (A-4)  into  equation  (A-l)  yields 

i  Bx  +  j  By  +  k  Bz  =  -(i  Mx  +  j  My  +  k  Mz)  /  R3 

+  3  (Mx  Rx  +  My  Ry  +  Mz  Rz)(  i  Rx  +  j  Ry  +  k  Rz)  /  R5.  (A-5) 

Collecting  terms  along  each  of  the  axes  yields 

iBx  =  -iMx/R3+3(MxRx  +  MyRy  +  MzRz)  iRx/R5  (A-6) 

j  By  =  -j  My  /  R3+  3  (Mx  Rx  +  My  Ry  +  Mz  R^  j  Ry/  R5  (A-7) 

k  Bz  =  -k  Mz  /  R3  +  3  (Mx  Rx  +  My  Ry  +  Mz  Rz)  k  Rz/  R5.  (A-8) 

Therefore,  the  magnitudes  of  the  east,  north,  and  upward  components  are, 
respectively, 

Bx  =  3  [Mx  x  +  My  y  +  Mz  (z  -  h)]  x/  R5-  Mx  /  R3  (A-9) 

By  =  3  [Mx  x  +  My  y  +  Mz  (z  -  h)]  y  /  R5  -  My  /  R3  (A-10) 

Bz  =  3  [Mx  x  +  My  y  +  Mz  (z  -  h)]  (z  -  h)/  R5-  Mz  /  R3.  (A-ll) 

The  anomaly  T  at  point  P  can  then  be  calculated  from 

T  =  (Bx  BEx  +  By  BEy  +  Bz  Bez)  /  Be  (y).  (A-12) 
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APPENDIX  B 


SEQUENCE  OF  STEREO  PAIR  REPRESENTATIONS  OF 
ANOMALY  FROM  A  POINT  FERROMAGNETIC  DIPOLE 


Thirty-six  different  (of  a  total  of  40)  horizontal  views  of  the  anomaly  from  a  point 
ferromagnetic  dipole  in  submarine-fixed  coordinates  illustrated  in  figure  4  are  shown  in  figure 
B-l.  The  angle  at  which  the  first  image  is  viewed  is  toward  the  east  at  095°  (magnetic)  and 
successive  images  are  viewed  at  10-degree  angular  increments,  proceeding  around  the 
compass  in  a  counterclockwise  manner.  Thus,  in  the  first  row,  the  viewing  angles  are  095°, 
085°,  075°,  065°,  055°,  045°,  035°,  025°,  015°,  and  005°.  The  image  at  the  end  of  each  row 
is  repeated  at  the  beginning  of  the  next  row.  Therefore,  in  the  second  row,  the  viewing  angles 
are  005°,  355°,  345°, . . .  275°.  The  ranges  of  angles  for  the  third  and  fourth  rows  are 
275° . . .  185°,  and  185° . . .  95°,  respectively.  The  images  are  arranged  such  that  they  can  be 
viewed  as  a  succession  of  stereo  pairs. 

The  submarine  was  assumed  to  be  on  a  north  magnetic  heading  in  the  Strait  of 
Hormuz  where  the  inclination  of  the  earth's  magnetic  field  is  40.3°  (downward).  The 
submarine  was  assumed  to  behave  as  a  point  ferromagnetic  dipole  having  a  resultant  moment 
of  4.39  x  107  y  ■  m3  (15.5  x  108  y  •  ft3)  inclined  downward  at  an  angle  of  27.2°.  The  axis  of 
highest  symmetry  of  the  anomaly  lies  in  the  plane  defined  by  the  earth's  field  vector  and  the 
magnetic  moment  vector  of  the  submarine  (in  this  case,  a  vertical  plane)  and  it  bisects  the 
angle  formed  by  these  two  vectors;  i.e.,  its  depression  angle  is  33.75°. 

The  toroid-shaped  portion  of  the  anomaly  has,  from  equation  (15),  an  outer  radius  of 
1140  m  (3741  ft)  and  its  surface  is  the  locus  of  points  for  which  the  anomaly  has  a  value  of 
-0.03  y.  Similarly,  from  equation  (14),  the  maximum  radial  extent  of  the  dumbbell-shaped 
portion  of  the  anomaly  is  1427  m  (4682  ft)  and  its  surface  is  the  locus  of  points  for  which  the 
anomaly  has  a  value  of  +0.03  y. 
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APPENDIX  C 


MAGNETIC  ANOMALY  ABOVE  AND  BELOW  THE  AIR-SEA 
INTERFACE  FROM  A  STATIC  HORIZONTAL  ELECTRIC  CURRENT 

ELEMENT  IN  SEAWATER 


Kraichman  (reference  (c))  provides  equations  that  enable  one  to  calculate  the  magnetic 
field  intensity  H  (in  amperes/meter)  in  each  of  two  semi-infinite  media  (e.g.,  atmosphere  and 
sea)  from  a  horizontal  electric  current  dipole  (HED)  in  the  lower  medium  positioned  on  the 
z-axis  at  an  elevation  h  relative  to  the  interface,  parallel  to  the  x-axis  and  positively  directed. 
(Note  that  h  requires  a  negative  sign  for  an  immersed  source.)  See  figure  C-l. 


Figure  C-l.  Geometry  and  equations  for  the  magnetic  field  intensity  above  and  within  a 
semi-infinite  conducting  medium  in  which  a  static  horizontal  electric  current 
dipole  is  immersed  (from  Kraichman,  1976). 
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Kraichman  defines  two  distance  parameters  Ri  and  R2  in  cylindrical  coordinates  as  the 
distances  to  a  field  point  P(p,  <Ji,  z)  from,  respectively,  the  current  element  in  the  lower 
medium  and  its  image  directly  above  it  in  the  upper  medium.  Figure  C-l  illustrates  the 
geometry  and  gives  Kraichman’s  equations  for  three  orthogonal  components  of  the  magnetic 
field  intensity  in  the  lower  medium  (medium  1)  and  the  upper  medium  (medium  2).  Note  that 
the  equations  for  the  vertical  component  Hz  are  identical  for  both  media  but  differences  in 
signs  occur  in  the  equations  for  the  other  two  components.  The  steady  electric  current 
moment  ps  is  the  product  of  the  steady  effective  current  Is  and  the  effective  length  L,  i.e., 
ps=  Is  L.  The  electrical  conductivities  of  medium  1  and  medium  2  are  ai  and  02,  respectively. 
If  the  earth’s  atmosphere  is  taken  as  medium  2,  its  conductivity  may  be  taken  as  zero;  this 
assumption  causes  all  of  the  factors  involving  at  and  02  to  vanish. 

In  the  absence  of  magnetic  materials  at  field  point  P,  one  can  convert  magnetic  field 
intensity  H  (in  A/m)  to  magnetic  induction  B  (in  nanotesla  or  gamma)  as  follows: 

1  A/m  =  4ji  x  10”7  T  =  400a;  nT  =  400ji  gamma.  (C-l) 

Thus,  the  equations  of  figure  C-l  can  be  recast  to  yield  the  magnetic  induction  at  field 


point  P  as  follows. 

Case  I.  z  <  0  (P  in  seawater) 

Bpi  =  -1 00  ps  sin  <j>  {(z-h)/Ri3  -  [(z+h)/R2  +1]  /p2}  (C-2) 

B*,  =  -100  ps  cos  <t>  {(z-h)/Ri3  +  (z+h)/R23  +  [(z+h)/R2  +  l]/p2}  (C-3) 

Bzi  =  100  ps  sin  <j>  p/Ri3  (C-4) 

Case  II.  z  >  0  (P  in  air) 

Bp2  =  -100  ps  sin  <j>  {(z-h)/R,3  +  [(z-h)/R,  -  1]  /p2}  (C-5) 

B*2  =  100  ps  cos  $  [(z-h)/Rj  -  l]/p2  (C-6) 

Bz2  =  100  ps  sin  <)»  p/R,3  (C-7) 


Case  III.  z  =  0  and  h  =  0  (Both  field  point  P  and  HED  at  air-water  interface) 


Bpi  =  Bp2  =  Bp  =  100  ps  sin  <J>  /  p2 

(C-8) 

B^i  =  B*2  =  B*  =  -100  ps  cos  4»  /  p2 

(C-9) 

Bzl  =  Bz2  =  Bz  =  100  ps  sin  <J>  /  p2 

(C-10) 
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For  calculating  the  magnetic  anomaly,  it  is  convenient  to  convert  to  Cartesian 
coordinates  and  then  to  compass  coordinates. 

Let  P(x,y,z)  =  field  point  at  which  anomaly  T  is  to  be  calculated 

x  =  distance  east  of  HED  (m) 

y  =  distance  north  of  HED  (m) 

z  =  distance  above  sea  level  (m) 

B,  =  easterly  component  of  magnetic  induction  from  HED  (y) 

By  =  northerly  component  of  magnetic  induction  from  HED  (y) 

Bz  =  vertically  upward  component  of  magnetic  induction  from  HED  (y) 

Be,  =  easterly  component  of  earth’s  magnetic  induction  (nT  or  y) 

Bey  =  northerly  component  of  earth’s  magnetic  induction  (nT  or  y) 

Bez  =  vertically  upward  component  of  earth’s  magnetic  induction  (nT  or  y) 

Be  =  magnitude  of  the  earth’s  magnetic  induction  vector  at  P  (nT  or  y) 
a  =  HED  compass  direction  (degrees  true) 

P  =  (x2  +  y2)1/2  (m) 

0  =  arctan  y/x  (degrees) 

<J>  =  a  +  0  -  90  (degrees) 

Be  =  (Be,2  +  BEy2  +  BEz2)1/2  (nT  or  y) 

For  either  Case  I,  Case  II,  or  Case  III 

Bx  =  Bp  cos  0  -  B„  sin  0  (y),  (C-l  1) 

By  =  Bp  sin  0  +  B*  cos  0  (y)  and  (C-12) 

T  =  (Bx  Bex  +  By  BEy  +  B,  Bez)  /  Be  (y)  (C-13) 
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APPENDIX  D 


MAGNETIC  ANOMALY  ABOVE  AND  BELOW  THE  AIR-SEA 
INTERFACE  FROM  AN  ALTERNATING  HORIZONTAL  ELECTRIC 
CURRENT  ELEMENT  IN  SEAWATER 


Kraichman  (reference  (c))  provides  equations  that  enable  one  to  calculate  the  magnetic 
field  intensity  H  (in  amperes/meter)  in  each  of  two  semi-infinite  media  (e.g.,  atmosphere  and 
sea)  from  a  horizontal  electric  current  dipole  (HED)  in  the  lower  medium  positioned  on  the 
z-axis  at  an  elevation  h  relative  to  the  interface,  parallel  to  the  x-axis  and  positively  directed. 
(Note  that  h  requires  a  negative  sign  for  an  immersed  source.)  See  figure  D-l  and  the 
drawing  of  figure  C-l.  The  distance  parameter  R  is  defined  in  cylindrical  coordinates  as  the 
distance  to  a  field  point  P(p,  <j>,  z)  from  a  point  on  the  air-sea  interface  directly  over  the 
alternating  HED. 


Subsurface  to  free-space  propagation  formulas  for  the  quasi-static  range, 

with  |y1R|»l  andR»  |h| 


Dipole 

tvpe 

H , 

<p 

»z 

HED 

p  sin  4>  eyi4  f„  3z2l 

2"y,  R3  l  R2J 

p  cos  <t>  eyi4 

2nyl  r3 

3p  sin  0  pey >*  f,  .  _  5z2l 

2  rry]  R5  L  >'1  R2J 

Subsurface  to  subsurface  propagation  formulas  for  the  quasi-static  range, 
when  \y0p\  «1  «|y1p|  with  p»|h|  and  |z| 

Dipole 

lype 

HED 

psin 

"ytp3 

pcos  *cy. (/>+*) 

2iry,p3 

3 P 

2 ny\pA 

Figure  D-l.  Equations  for  (a)  subsurface-to-free-space  propagation  and  (b)  subsurface-to- 
subsurfacc  propagation  for  the  quasi-static  range  (from  Kraichman,  1976). 


The  symbols  used  in  figure  D-l  are  the  same  as  those  used  in  Appendix  C  except 
p  =  p,  =  the  alternating  electric  current  moment 

R  =  (p2  +  z2)m  =  distance  to  a  field  point  P(p,  <(>,  z)  from  a  point  on  the  air-sea 
interface  directly  over  the  alternating  HED 
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Yo  =  propagation  constant  for  free  space  (air) 

Yi  =  propagation  constant  for  a  homogeneous  conducting  half-space  (seawater). 
Kraichman  gives  the  propagation  constant  y  through  the  expression 

Y2  =  -co2p  e  +j  cop  a  (D-l) 

in  which  to  =  the  angular  frequency  of  the  source  =  2  ji  f 

p  =  the  magnetic  permeability  constant  of  the  medium 
e  =  the  electrical  permittivity  constant  of  the  medium 

j  =  (-1)1* 

a  =  the  electrical  conductivity  of  the  medium. 

For  air  or  seawater,  p  =  po  =  4  n  x  10-7  T  •  m/A.  In  air,  e  =  £0  =  8.85  x  lO-12  C?/N  •  m2; 
in  seawater,  £  =78  eo-  For  air,  ct  =  02  =  0  and,  for  seawater,  a  =  Oj  =  4  S/m.  In  the  following 
discussion  of  the  conditions  under  which  the  equations  listed  in  figure  D-l  are  valid  it  is 
assumed  that  the  frequencies  of  interest  are  less  than  100  Hz. 

Field  Point  Pin  Air 

Condition  1.  |yi  R|  »  1 

For  frequencies  of  interest,  it  can  be  shown  that  the  first  term  on  the  right  side  of 
equation  (D-l)  is  many  orders  of  magnitude  smaller  than  the  second  term  and,  therefore, 
negligible  in  comparison.  Thus,  the  modulus 

|Yi|  =  (cu  p  a) 1/2  -  (2  ji  f  po  d)m  -  5.620  x  10"3  f1/2  or  (D-2) 

5.620  x  lO'3  fw  R  »  1.  (D-3) 

It  follows  that  the  range  R  must  therefore  satisfy  the  condition 

R  »  177.9/f,/2.  (D-4) 

If  the  symbol  “»”  is  interpreted  as  meaning  “ten  times  greater”  then  the  range  R  must  exceed 
1779/f1/z.  Thus,  for  f  =  1  Hz,  R  must  exceed  1779  m  and,  for  f  =  4  Hz,  R  must  exceed  890  m. 
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Condition  2.  R  »  |h| 


If,  for  example,  the  depth  of  the  HED  is  100  m,  the  equations  are  valid  for  ranges  R 
greater  than  1000  m,  if  the  above  interpretation  of  “»”  is  held.  Because  R  is  reckoned  from 
the  sea  surface  (rather  than  from  the  HED),  it  is  actually  possible  for  R  to  be  less  than  the 
depth  |h|  of  the  HED.  Thus,  care  must  be  exercised  to  insure  that,  when  calculations  are 
performed  for  P  nearly  or  directly  over  the  HED,  the  altitude  of  P  is  much  larger  than  the 
depth  of  the  HED. 

Field  Point  P  in  Seawater 

Condition  1.  |y0  p|  «  1 

Because  the  conductivity  of  air  a  =  o2  =  0,  the  second  term  of  equation  (D-l)  vanishes 
when  the  modulus  |y0|  of  the  propagation  constant  in  air  is  calculated.  Thus, 

|y0|  =  co(poe0)1/2  -  2  n  f  (po  eo)m  -  2  n  f  /c  =  2.096  x  10~8  f  (nT1)  (D-5) 

and  therefore  the  condition 


|yo  p|  =  |2.096  x  10'8  f  p|  «  1  (D-6) 

is  satisfied  easily  for  the  frequencies  f  and  ranges  p  of  interest. 

Condition  2.  |yi  p|  »  1 

For  frequencies  of  interest,  it  can  be  shown  that  the  first  term  on  the  right  side  of 
equation  (D-l)  is  many  orders  of  magnitude  smaller  than  the  second  term  and,  therefore, 
negligible  in  comparison.  Thus,  the  modulus 

|Yi|  =  (w  p  o)m  -  (2  rc  f  po  o)m  =  5.620  x  10'3  f1/2  (m’1)  or  (D-7) 

5.620  x  10"3  tm  p  »  1.  (D-8) 

It  follows  that  the  horizontal  component  of  range  p  must  therefore  satisfy  the  condition 

p » 177.9/f1/2  (m).  (D-9) 

If  the  symbol  “»”  is  interpreted  as  meaning  “ten  times  greater”  then  the  range  p  must  exceed 
1779/fi/2  Thus,  for  f  =  1  Hz,  p  must  exceed  1779  m  and,  for  f  =  4  Hz,  p  must  exceed  890  m. 

Condition  3.  p  »  |h| 

If,  for  example,  the  depth  of  the  HED  is  100  m,  the  equations  are  valid  for  ranges  p 
greater  than  1000  m,  if  the  above  interpretation  of  “»”  is  held. 
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Condition  4.  p  »  |z 


If  the  above  interpretation  of  “»”  is  held,  the  equations  are  valid  only  if  the  horizontal 
component  of  range  p  of  the  underwater  field  point  P  exceeds  ten  times  the  depth  of  P;  i.e., 
only  if  the  angle  subtended  at  a  point  on  the  surface  directly  above  the  HED  between  the  sea 
surface  and  point  P  is  less  than  5.7°  (i.e.,  arctan  0.1).  Thus  the  equations  of  figure  D-l  cannot 
be  validly  applied  to  the  full  hemisphere  below  the  HED. 

Subject  to  the  above  six  conditions,  the  equations  of  figure  D-l  can  be  recast  to  yield 
the  magnetic  induction  at  field  point  P  in  a  manner  similar  to  that  used  in  Appendix  C  as 


follows. 

Case  I.  z  <  0  (P  in  seawater) 

Bpi  =  [400  pa  sin  <j>  /  (yi  p3)]  {expfyi  (h  +  z)]}  (D-10) 

B*i  =  [-200  pa  cos  <J>  /  (yi  p3)]  {exp[yi  (h  +  z)]}  (D-l  1) 

Bzi  =  [600  pa  sin  <f»  /  (yi2  p4)]  (exp[Yi  (h  +  z)]}  (D-12) 

Case  II.  z  >  0  (P  in  air) 

Bp2  =  [200  Pa  sin  <j>  /  (Yr  R3)]  [exp(Yi  h)]  [2  -  (3  z2/  R2)]  (D-13) 

B<,2  =  [-200  p.  cos  <j>  /  (yi  R3)]  [exp(Yi  h)]  (D-14) 

Bz2  =  [600  p,  sin  <j>  /  (yi2  R5)]  [p  exp(Yi  h)]  [1  +  Yi  z  -  (5  z 2 1 R2)]  (D-15) 

Case  III.  z  =  0  and  h  =  0  (Both  field  point  P  and  HED  at  air-water  interface) 

Bp  =  400  pa  sin  <j>  /  (yi  p3)  (D-16) 

B*  =  -200  p,  cos  <j>  /  (yi  p3)  (D-17) 

Bz  =  600  p.  sin  <j>  /  (yi2  p4)  (D-18) 


Equations  (C-ll),  (C-12),  and  (C-13)  can  be  used  with  the  preceding  nine  equations 
to  determine  the  magnetic  anomaly  produced  in  the  earth’s  field  in  each  of  the  above  three 
cases. 
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